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' We present and discuss some ideas concerning an "average-pair-density functional theory", in 

which the ground-state energy of a many-electron system is rewritten as a functional of the spheri- 
cally and system-averaged pair density /(r^). These ideas are further clarified with simple physical 
examples. We then show that the proposed formalism can be combined with density functional 
theory to build system-adapted correlation energy functionals. A simple approximation for the un- 
known effective electron-electron interaction that enters in this combined approach is described, and 
results for the He series and for the uniform electron gas are briefly reviewed. 
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Density Functional Theory (DFT) is nowadays the 
most widely used method for the calculation of electronic 
structure in both solid-state physics and quantum chem- 
istry [H Q ■ ^ e accurac y °f the results coming from 
a DFT calculation is limited by the approximate nature 
of the exchange and correlation energy E xc [n] and the 
associated potential, its functional derivative. While, at 
present, many algorithms for a very accurate or even ex- 
act exchange potential v x (r) and energy E x [n] are avail- 
able |ij.|iJ]- better approximations for calculating accurate 
correlation energies E c [n] are still needed. 

In a recent paper [5j, we proposed an alternative ap- 
proach to build the DFT correlation energy E c [n\. The 
method consists in solving simple radial equations to gen- 
erate the spherically and system-averaged pair density 
/(T12) along the so-called adiabatic connection. Besides 
its practical use for DFT, we realized that this method 
has many aspects that deserve to be better investigated, 
including the possibility for an alternative theory com- 
pletely based on /(Y12). In this work, we sketch the basic 
ideas for this alternative theory, we further discuss them 
with simple physical examples, and we deal with some of 
the aspects that went overlooked in Ref. Jjj. The scope 
of this paper is to lay the rigorous foundations for an 
approach that will be further developed towards the con- 
struction of a self-consistent scheme combining the radial 
equations for f{r\2) with the Kohn-Sham equations. 

The paper is organized as follows. After defining the 
notation, we develop in Sec. II I II the formalism corre- 
sponding to a theory based on the spherically and system 
averaged pair density, reviewing at the same time the cor- 
responding concepts for DFT. We hope that this parallel 
treatment will make it easier for the reader to familiar- 
ize with the new concepts. We then give, in Sec. IIVI 
some simple physical examples. Section [V] is devoted to 
explain how the ideas of Sec. II I II can be used to build 
correlation energy functionals for DFT. A simple, physi- 
cally motivated, approximation for the unknown effective 
eletron-electron interaction that appears in our formal- 
ism is discussed in Sec. I VII where applications to two- 
electron atoms and to the uniform electron gas are also 
briefly reviewed. The last section is devoted to conclu- 



II. DEFINITIONS 

We start from the standard Af-electron hamiltonian 
(in Hartree atomic units, h = m = a = e = 1, used 
throughout) 



H = 
T = 

V,, = 



T + V ee + V ne , 
1 N 

i=l 
1 N 1 

N 



(1) 

(2) 
(3) 
(4) 



where v ne is the external potential due to nuclei. Given 
the ground-state wavefunction of H, we consider two 
reduced quantities that fully determine, respectively, the 
expectation values (^\V ne \^) and (^\V ee \^), i.e., the 
usual one-electron density 

n(r) = N J \y{r,r 2 ,...,r N )\ 2 dr 2 ...dr N , (5) 



and the spherically and system-averaged pair density 
(APD), which is obtained as an integral of | \[/ 1 2 over all 
variables but r i2 = |i"i — r 2 |, 



/(ri 2 ) 



N(N — 1) 



E 



|*(r 12 ,R, r 3 ,...,r N )\ 2 ^^dRdr 3 ...dr N , 



(6) 



where ri 2 = ri — r 2 , and R = |(ri + r 2 ). The func- 
tion JJrrz) is also known in chemistry as intracule den- 
sity (6L F7l Isl I9L fiol fill TT^I . and, when multiplied by the 
volume element 4irr 2 2 dr 12, is proportional to the proba- 
bility distribution for the electron-electron distance. We 
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then have 



B. APDFT — The system-dependent functional 



(*|Ke|*> = / n(r)v ne (v)dr (7) 

_ r /M^^. (8) 



ri2 



-dr i2 



i: ri2 



In the following text we will also deal with modified sys- 
tems in which the external potential and/or the electron- 
electron interaction is changed. Thus, the notation V ee 
and V ne will be used to characterize the physical system, 
while the modified systems will be defined by W and V, 
with 



1 N 



N 



v = 5>(r0, 



(9) 



(10) 



»=i 



where the pairwise interaction w always depends only on 



III. FORMALISM 

In this section we present a "APD-functional theory" 
(APDFT) based on the function /(m) highlighting, step 
by step, the analogies in reasoning with the derivation of 
standard DFT. 



A. DFT — The universal functional 

In standard DFT one defines a universal functional of 
the one-electron density n as resulting from a constrained 
search over all the antisymmetric wavefunctions \& that 
yield n [13 

F[n;V ee) T\ = wm(*\T+V ee \*), (11) 
or, more completely, as a Legendre transform |14| 
F[n;V ee ,T] = sup { nrin<*|T + V ee + V\V) 

i(r>(r)<2r}. (12) 

In both Eqs. (|ll|l and 1|12|1 . the dependence on the 
electron-electron interaction (and on the kinetic energy 
operator T) has been explictly shown in the functional. 
The universality of the functional F stems exactly from 
the fact that the e-e interaction is always 1/r (and that 
T is always the same). 

The ground-state energy Eq of the system can then be 
obtained by minimizing the energy with respect to n, 



E =min^F[n;V ee ,T] + J 



n(r)v ne (r)dr 



(13) 



Similarly, we can define a system-dependent functional 
(i.e., a functional depending on the external potential 
V ne , and thus on the specific system) of the APD f(r 12) 
as 



G[f;V ne ,T\ = mm(*\T + V ne \V), (14) 



or better as 



G[f;V ne ,T] = Bup{min<*|T + W + V Be |tf) 

f{ri 2 )w(ri 2 )dr 12 \. (15) 



The ground-state energy can be obtained by a minimiza- 
tion with respect to / 



E Q = min (g[/; V ne ,T] + f ^1 dr 12 
f { J r 12 



(16) 



Evidently, with respect to DFT, the functional G has the 
disadvantage of being not universal: in DFT, an approx- 
imation for F should be in principle valid for all systems. 
However, the crucial point for applications is understand- 
ing how difficult is to build a reasonable approximation 
for G[f;V ne ,T], given a certain V ne . In the particular 
combination of DFT and APDFT that we propose in 
Sec. |3 the lack of universality of G is not an issue. 

Another issue concerns the A-representability condi- 
tions on /(r 12 ), i.e., which contraints must a given / 
satisfy to guarantee that it comes from the contrac- 
tion of a A-electron wavefunction , 3>. This is a prob- 
lem shared with other generalizations of DFT, like the 
pair-density functional theory [H H [H H El The 
iV-representability conditions on f(r 12) are evidently re- 
lated to those on the pair density, and we thus might ex- 
pect that they are not a trivial matter [Isj . The definition 
of the functional G of Eq. 115|) formally overcomes this 
problem by giving a divergent (+00) answer for any non 
iV-representable / (when / is not A-representable the 
right-hand side of Eq. (|15|l is not bounded from above) , 
but this is not a practical solution when coming to appli- 
cations. As we shall see, in the particular use of APDFT 
presented in Sec. this issue is not particularly crucial 
from a practical point of view, since we use APDFT to 
build correlation functionals for DFT 



C. DFT — Adiabatic connection 

In density functional theory, one usually defines a set 
of hamiltonians depending on a parameter A [2(], El , 



H = T + W + V 



(17) 



having all the same one-electron density, equal to the one 
of the physical system 

n A (r)=n(r) VA. (18) 
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If W x=0 = and W XphyB = V ee , one switches continu- 
ously from a noninteracting system to the physical sys- 
tem, while keeping the density fixed by means of a suit- 
able external potential V x . Obviously, the APD f(r 12) 
changes with A. By the Hellmann-Feynmann theorem, 



dX 



dW x dV x 

'^T J — 



dX 



A ,dw x (r 12 ) 



(19) 



so that by directly integrating Eq. (|19|l , and by combining 
it with Eq. JT^, one obtains 



F[n; V ee ,T] = T s [n] - 



dX / dr 12 f x (r 12 



dw > 



J12 



dX 



where T s [n] is the kinetic energy of a noninteracting sys- 
tem of N electrons with density n(r). 

More generally, one can be interested in using as a 
starting point a system of partially interacting electrons, 
corresponding to a particular value of the coupling A (say, 
A = /1) between and A p h ys - In this case, if is the 
wavefunction of the system with partial interaction 
(and external potential V 1 - 1 ) we have 



F[n;V ee ,T] 



+ 



dX 



dr 12 f\r 12 ) d -^21) 



Usually, the adiabatic connection is performed along a 
linear path by setting W x = XV ee (thus A p h ys = 1), 
but some nonlinear choices can be more convenient when 
dealing with approximations. 



so that 



G[f;V ne ,T]=T { [f] + 



£phys 



d£ / drn^r) 



dv^ (r) 



(25) 



where Tf[f] is the kinetic energy of a system of N free 
fermions (zero external potential) having the same f{r\ 2 ) 
of the physical system. A simple example of such adi- 
abatic connection is given in Sec. IIV Al As we shall 
see, given a confined physical system, the correspond- 
ing w^ =Q (r 12) must be partially attractive (in order to 
create a bound cluster of fermions). This could in princi- 
ple lead to "exotic" ground states for some of the systems 
corresponding to £ < £ p h ys - This issue is not considered 
in this paper, and will be investigated in future work. 

Similarly to the DFT case, it could be convenient to 
choose as starting point a system with an external po- 
tential corresponding to some coupling constant £ (say, 
£ = a) between and £ p h ys - If ^ a is the ground-state 
wavefunction of the system with external potential V a 
(and e-e interaction W a ) we have 

G[f;V ne ,T] = (y a \T + V a \V a ) + 

dtj drn«(r)^ii. (26) 



E. DFT — Kohn-Sham equations 

One-particle equations in DFT can be obtained by 
defining a set of orthogonal orbitals <Pi(r) with occupa- 
tion number Vi that minimize X^i^^il — \^ 2 \ L Pi) an d 
yield the density of the physical system, J^i v i\ l fii{ v )\ 2 = 
n(r). This gives 



D. APDFT - Adiabatic connection 



We can also define a set of hamiltonians 



= T + W^ + 



(22) 



in which the function f{r\ 2 ) is kept fixed, equal to the 
one of the physical system, 



fHru) = /(ria) V£, 



(23) 



If V^ =0 = and V^ phys = V ne , we are switching con- 
tinuously from a system of N free electrons interacting 
with a modified potential w^ =0 (ri 2 ), to the physical sys- 
tem. That is, f(ri 2 ) is kept fixed as £ varies by means 
of a suitable electron-electron interaction while the 
one-electron density n(r) changes with £. Again, by the 
Hellmann-Feynmann theorem, we find 



' ,dw^(r 12 ) 
/(ri 2 ) ^ dr 12 



rz> (r) 



dr, (24) 



[-iy 2 + Vl (r)] ipi(r) = e lVl (r) 

5>|^(r)| 2 = n(r), (27) 

i 

where the potential t>i(r) is the Lagrange parameter for 
the density. To fully specify these equations one needs 
a rule for the occupation vi of the orbitals. The Kohn- 
Sham choice corresponds to occupy the orbitals in the 
same way as for a Slater determinant. This determinant 
is the wavefunction of a system of N non- interacting elec- 
trons constrained to have the same one-electron density 
of the physical system, and leads to the identification 



T s [n] — min 



X 



(28) 



with the same T s [n] of Eq. (121 H . The ground-state energy 
of the physical system is then obtained via the Hartree- 
cxchange-correlation functional E-R XC [n], defined as the 
difference F[n; V ee ,T] — T s [n]. This also implies that, in 
Eqs. Vi(r) = v KS (r) = v ne (r) + SE Hxc [n]/8n(r). 

Usually, the Kohn-Sham equations are derived starting 
from the noninteracting system with density n(r), rather 
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than from a constrained minimization of Vi (ifi 



±V 5 



This different way of proceeding allows us to 



keep the analogy with what we will do in the next sub- 
section for APDFT. 



F. APDFT - effective equations 

Since the e-e interaction is spherically symmetric, the 
relevant APD that determines (^|T4 e |^) is a unidimen- 
sional quantity. To obtain simple "two-particle" equa- 
tions for f{r 12) we start from the kinetic energy operator 
for the scalar relative coordinate r\ 2 = Ir2 — ri|, 



ri2 dr 



£{£ + !) 



(29) 



12 



12 



and we perform a minimization of X^^i^il — ^'r^lV'i) 
with respect to some orthogonal "effective" geminals 
V'i( r i2) constrained to yield / of the physical system, 
Y^i^Mi? = fi leading to 



[-Vj? 12 + w cS (ri 2 )}ipi(r 12 ) = eiipi(r 12 ) 



(30) 



- T g [f}. This also 
HP, w cS (r 12 ) = 



the kinetic and external-potential functional defined in 
Ref. 0, F KE [f;V ne ] = G[f;V ne ,T] - 
leads to the identification, in Eqs 
l/r 12 + 5F KE [f;V ne }/Sf(r 12 ). 

An important issue to be addressed concerning the ef- 
fective equations (|30fl is whether a given physical (and 
thus iV-representable) f(r 12) is also representable by the 
simple "effective-geminal" decomposition of Eqs. J3DJ. 
This question is similar to the one arising in DFT: is 
a physical density always non-interacting representable? 
In view of the more complex nature of f{r\ 2 ) with re- 
spect to n(r) we might expect that this problem is much 
more difficult to face in APDFT than in DFT. It seems 
reasonable that at least the short-range part of a physical 
f( r i2) is representable by Eqs. I|3l)fl . while the long-range 
tail of / of an extended system could be problematic [26| . 



IV. SIMPLE PHYSICAL EXAMPLES: 
TWO-ELECTRON SYSTEMS 

In this section we give some examples for two-electron 
systems, in order to gain physical insight with some of 
the ideas just introduced. 



The interaction w e s(ri 2 ) is the Lagrange parameter for 
/. Again, to fully specify these equations we need a rule 
for the occupancy $i of the effective geminals. For spin 
compensated systems, we can choose to apply a rule that 
resembles to a Slater determinant: occupancy 1 for even 
£ (singlet symmetry), occupancy 3 for odd £ (triplet sym- 
metry), up to N(N — l)/2 occupied geminals. This rule 
has been applied to solve the effective equations H3U|) 
in the uniform electron gas, with rather accurate re- 
sults [23l IziL |25| . It is however important to point out 
that when we apply this occupancy rule to Eqs. I|3U|) 

1. there is no Slater determinant that can be asso- 
ciated with our effective geminals: the tpi are con- 
strained to give the exact f that cannot be obtained 
from a noninteracting wavcfunction (for example, 
any Slater determinant violates the cusp condition 
satisfied by the exact /); 

2. more generally, there is no wavefunction (and so no 
physical system) that can be built from our effective 
geminals. 



This last point implies that, if we define 

T g [f]= min Y>i|-V? 12 |Vi 



(31) 



(with the determinant-like occupancy), we have in gen- 
eral 



T g \f\*Tt\f\, 



(32) 



where Tf[f] was defined in Eq. The total en- 

ergy of the physical system can then be recovered via 



A. A picture from harmonic forces 

A very simple picture of the whole adiabatic connection 
path in both DFT and APDFT can be gained by looking 
at an analytic-solvable model, i.e., a two-electron hamil- 
tonian with only harmonic forces (harmonic electron- 
nucleus attraction, and harmonic e-e repulsion too) [TTj . 



H{K,k) = -\Vl-\Vl + \Kr 2 l + \Krl-\k\v 1 -v 2 \\ 

(33) 

where K > (attractive nucleus-electron potential), 
k > (repulsive c-c interaction), and K > 2k, to have a 
bound system. For this hamiltonian, 



n(r) 



f(ri2) 



20 



3/2 



T 3/2 
y 3/2 
TV? 



-fir 2 - 



~T r 12 



a 2VKW^W) 
^K-2k + VK 



7 = -VK -2k. (35) 



1. DFT 



If our "physical" system corresponds to some K = K n 



and k 



and we want to switch off the e-e interaction 



(by setting, e.g., k = Xk ee ) while keeping the density 
fixed, we will simply have to change K ne into K(\) such 
that in Eq. l|3^|l does not change. The function K(X) is 
shown in the upper panel of Fig. ^ for the case K ne = 3, 
k ee = 1. We see that, as A — * 0, K(X) decreases, because 
a smaller attraction is needed to keep the electrons in 
the density when there is no e-e repulsion. Of course, 
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f(r±2) changes with A, as shown in in the lower panel of 
Pig. Q] as A decreases, the "on-top" value / A (r 12 = 0) 
gets larger. This reflects the fact that when there is no 
e-e repulsion, it is more likely to find the two electrons 
close to each other. 



2. APDFT 

If, instead, we switch off the external potential, (e.g., 
by setting K — £K ne ) while keeping f(r 12) fixed, we will 
have to change the e-e interaction in order to keep 7 of 
Eq. i|35|) constant, 



m = w~i)K ne + k e 



(36) 



The one-electron density along this adiabatic connection 
is 



-0(0 r 2 



7T 3 /2 

2 \J £,K ne (K ne 1k ee ) 



y/K ne - 2k ee + vW„ 



(37) 



Thus n^(r) is a gaussian that, as £ — ► 0, becomes more 
and more spread, as shown in the lower panel of Fig. 
When the external potential goes to zero, the system 
becomes translationally invariant and the wavefunction 
for the center-of-mass degree of freedom is simply a plane 
wave. Correspondingly, the electron-electron interaction 
changes with Eq. ffify : we see from the upper panel of 
Fig-Glthat, as £ gets smaller, £;(£) becomes smaller (less 
repulsive), and then it changes sign at some < £ < 
1, becoming an attractive interaction. For a confined 
system, when the external potential approaches zero, an 
attractive e-e interaction is needed in order to keep f(r 12) 
fixed. Moreover, in the very special case of harmonic 
forces there is a value £* = l—2k e e/K ne £ (0, 1) for which 
the e-e effective interaction is zero everywhere, tu* = 0. 

Of course, when the e-e physical interaction is the 
Coulomb potential l/ri2 this cannot happen: a system 
with the same f(r 12) of the physical system cannot have 
w = everywhere. This can be simply understood by 
thinking that there is no external potential that can force 
the system to have the correct cusp 27] at rn = 0. In 
fact, along any adiabatic connection that keeps f(r 12) 
fixed, equal to the one of a system with Coulombic e-e 
interaction, w^(r 12) will always behave as l/ri2 in the 
limit of small ryi- 



B. He atom 

Consider now the two electrons of a He atom. They 
feel the attraction of the nucleus, — 2/r, and they repel 
each other with potential l/ri 2 . Given the exact (or a 
very accurate |28|) ground-state wavefunction 5', we can 
calculate the "exact" density n(r) and the "exact" /(ri 2 ). 




'12 



FIG. 1: Adiabatic connection in DFT for the simple harmonic 
hamiltonian of Eq. 1331 . The e-e interaction is multiplied by 
a parameter A, and the density is kept fixed by a suitable 
external potential (upper panel). The APD f(rrz) changes 
with A as shown in the lower panel. 




FIG. 2: Adiabatic connection in APDFT for the simple har- 
monic hamiltonian of Eq. 1331 . The external potential is mul- 
tiplied by a parameter £, and the function /(ri2) is kept fixed 
by a suitable electron-electron interaction (upper panel). The 
density n(r) changes with £ as shown in the lower panel, and 
as £ — > becomes completely delocalized. 
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Now, we can consider the case in which W — and n(r) 
is kept fixed (the KS system in DFT, Sec. llll (Jjl . and the 
one in which V = and /(V12) is kept fixed (APDFT, 
SecHITDl): 



1. DFT 

We construct a system which has the same density n(r) 
of the physical one and no electron-electron interaction. 
This is the KS system, in which the two electrons do not 
interact (w = 0) and feel an external potential v(r) less 
attractive than — 2/r, as in the case of the harmonic po- 
tential of Fig. [I] The APD / A=0 (ri 2 ) of this system will 
be different from the physical one, as shown in Fig. [3] We 
see that the change in the function / when we switch from 
the physical system to the KS one is quantitatively simi- 
lar to the one of Fig. i.e., at A = the "on-top" value is 
higher than the physical one. In the case of Coulomb e-e 
interaction the physical f{r\2) has a cusp, f (0) = /(0), 
due to the short-range divergence of l/r\2 |27|. 



2. APDFT 

In APDFT, we can construct a system which has the 
same f{r 12) of the physical one, and zero external poten- 
tial (V = 0). This is a system of two bounded fermions 
interacting with the effective potential w^ =Q (r 12) of Fig.0] 
(calculated from an accurate [HEl] /). As in the case of 
harmonic forces, the density of this system is completely 
delocalized, because the wavefunction for the center-of- 
mass degree of freedom is a plane wave. We can imag- 
ine that along the linear adiabatic connection, v^(r) = 
— 2£/r, the corresponding w^(r 12) changes smoothly be- 
tween l/ri2 (at £ = £ p hys = 1) and the potential 
w^ =0 (r 12) of Fig. 0] As anticipated, we se that at £ = 
the effective e-e interaction has an attractive part, which 
is necessary to have the same / of a physical confined sys- 
tem. However, as already pointed out, for small ri2, the 
e-e interaction must always behave as l/ri2, to produce 
the exact cusp in /(ri 2 ) |27| . 



V. FROM APDFT TO CORRELATION 
ENERGY FUNCTIONALS FOR DFT 

In Sec. lIIII we have underlined the similarity of the roles 
played by n(r) and f(r 12) from a mathematical point of 
view: the former completely determines (\E , |V / „e|5'), and 
the latter (5 , |V r ee |^ / ). However, while the KS system is a 
substantial simplification of the many-electron problem 
(yielding to single particle equations), the auxiliary sys- 
tem of Sec. IIII Dl (with zero external potential) is still a 
complicated many-body object, consisting of N fermions 
interacting with a partially attractive potential. The ra- 
dial equations of Sec. IIII Fl are a great simplification of 




0.5 1 1.5 2 2.5 

r 12 



FIG. 3: The function f x (r\2) at the two ends of the adiabatic 
connection in DFT for the He atom. For A = 0, we have a 
system of two noninteracting (W = 0) electrons constrained 
by an external potential to yield the same density of the phys- 
ical system. For A = A p h ys we have the physical system with 
full interaction \jr\i and external potential — 2/r (from the 
wavefunction of Ref. |2j|; see also Ref. |f|). 
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FIG. 4: The electron-electron interaction vrirvi) at the two 
ends of the adiabatic connection in APDFT for the He atom. 
For £ = we have a system of two free fermions (zero ex- 
ternal potential) interacting with MT =0 (ri2). This system has 
the same /(ri2) of the physical system, but a completely de- 
localized one-electron density. For £ = £ p h ys we have the 
physical system, with e-e interaction l/r-12 and external po- 
tential —2/r. (The potential w at £ = has been calculated 
from the accurate wavefunction of Ref. |2fi|: see Ref. [f| for 
more details.) 



the problem, but we might expect that building approx- 
imations for the whole functional Fke[/; Vne] could be 
not easy. 

Our basic idea, instead, is to use APDFT to build what 
is missing in DFT, i.e., to build / A (ri 2 ) along the adia- 
batic connection in DFT. As said in the Introduction, we 
insert our approach in the framework of exact-exchange 
DFT 0, 0, 123. |3(j in which only the correlation energy 
functional needs to be approximated. The ground-state 
energy of the physical system is given by 

E = T s [n] + J n(r)v ne (r)dr + E H [n] + E x [n] + E c [n], 

(38) 

where -EhM is the usual Hartree term, -E x [n] is the ex- 
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change energy, obtained by putting the Kohn-Sham or- 
bitals in the Hartree-Fock expression for exchange, and 
the correlation energy, £7 c [n], is unkown. Equation 1(38(1 
can be also rewritten as 

E = ($ KS |T + V ee + V ne \$ KS ) + E c [n], (39) 

where $ks is the Slater determinant of Kohn-Sham or- 
bitals, i.e., the wavefunction of N noninteracting elec- 
trons constrained to yield the same n(r) of the physical 
system. Combining Eq. ((20(1 with Eq. 1)39(1 . we see that 
the wanted correlation energy is given by 

E c [n] = dX £ dr 12 4,r? 2 ^(r,,)^ , 

(40) 

where 

tiiru) = f\r 12 ) - / A =°(r 12 ) = f\r 12 ) - f KS (r 12 ). 

( 41 ) 

Thus, in order to get the KS correlation energy, we should 
compute f x (ri 2 ) for each hamiltonian H x of Sec. IIII CI 
Our approach consists in solving the simple radial equa- 
tions of Sec. IIII Fl for each H A along the adiabatic con- 
nection in DFT. This is not particularly expensive: we 
are dealing with unidimensional equations, and, if the de- 
pendence of w x on A is smooth, we will only need few A 
values (~ 5 — 30) between and A p h ys - With this partic- 
ular combination of APDFT and DFT, we do not need 
to approximate the whole functional Fke along the DFT 
adiabatic connection, but only its functional derivative, 
i.e., the effective interaction w c g(ri 2 ) which appears in 
Eqs. 1(30(1 . since the remaining information is provided 
by DFT. As we shall see, simple physical arguments can 
be used to build reasonable approximations for w e s{r\ 2 ) 
at each coupling strength A. 

Since the effective equations yielding f x (ri 2 ) must be 
solved for each system, we speak of system-adapted cor- 
relation energy density functionals. 

VI. BUILDING APPROXIMATIONS 

In Ref. we proposed and successfully tested a sim- 
ple approximation for building w c ^{ri 2 ) along the DFT 
adiabatic connection for two-electron atoms. This ap- 
proximation starts from w^{ri 2 ), the effective e-e inter- 
action that gives /ks( r i2) when inserted in Eqs. I|30|l . In 
the special case of two-electron systems, (ri 2 ) is di- 
rectly available in a KS calculation. For systems with 

more than two electrons, w^g (Vi 2) could be calculated, 
e.g., with the methods of Refs. tHHl^. Then, the idea 
is to build an approximation for a correlation potential, 

to be added to w^g (^12), which describes the change in 
/ when the e-e interaction is turned on, from zero to 
w x (r \ 2 ). To do this, we defined an average density n, 

n=l /drn(r) 2 , (42) 



- r, ? Wet' 


He . 

approx. 






,--""(0) w eff ' exact ' 
w eff 



0.5 1 1.5 2 

r 12 



FIG. 5: Construction of an approximation for the effective 
potential that generates the APD f(r\2) of the He atom: w^g 
is the part of the potential that generates the APD of the 
Kohn-Sham system. The "exact" [3, 28] correlation potential 
Wefr and our approximation of Eq. 1 lot are shown, together 
with the Coulomb repulsion l/ri2. 

and, correspondingly, an average radius r s , 

r s = (fn)- 1/3 . (43) 
We then built a correlation potential w^,g(ri 2 ), as 

w C eg(r 12 )=w x (r 12 )- I nw x (\r- r 12 \)dr. (44) 

J\r\<r B 

The idea behind Eq. ((44(1 is the following: the e-e in- 
teraction u> A is screened by a sphere of radius r s and 
of positive uniform charge of density n that attracts the 
electrons with the same interaction w A . The average den- 
sity n of Eq. ((42(1 (and thus the average r s ) is kept fixed 
to mimic the fact that the one-electron density does not 
change along the adiabatic connection. 

In order to gain insight with our construction, let us 
consider the case of the physical system, w A=A P h y= = 
l/ r i2> for which Eq. 1(441) corresponds to 

< ff (ri2) = ( — + % £-) (f s r 12 ) , (45) 
\ri 2 27* 2r s J 

where 9{x) is the Heaviside step function. In Fig. [3] we 
reported, for the He atom, the potential w^g which gen- 
erates /ks, together with the "exact" correlation poten- 
tial W^g, and the approximation of Eq. ((45(1 . We see that 

the potential w^g is a confining potential for the vari- 
able r±2'. our idea is to include in this term, available 
from DFT, the contribution to f{ri 2 ) coming from the 
particular external potential of the system and from the 
fermionic structure of the wavefunction. The remaining 
part to be approximated, the correlation potential w^g, 
must include the effect of the c-e repulsion while keeping 
the density fixed, i.e., it must be essentially a screened 
Coulomb interaction. We see from Fig. that the simple 
approximation of Eq. 1(45(1 is reasonable, i.e., the screen- 
ing length is well approximated by r s of Eqs. ((42|) - ((43(l . 



8 



For comparison, the full Coulomb repulsion 1 /ri 2 is also 
shown. Notice that in the special case of two-electron sys- 
tems we have T g [f] — Tf[f], so that the potential w^ =0 

of Fig. ^corresponds, in Fig. [SJ to the sum of and 
the "exact" w^ s . 

In Ref. we inserted the potential ucg (ri2) + 
w^g (7*12) into Eqs. (|30(l . and solved them for several two- 
electron atoms. Our results can be summarized as fol- 
lows: (i) at A = Aphys (i.e., for w x (r 12) = 1/7*12) we ob- 
tained APD f(r 12) in close agreement with those coming 
from accurate variational wavefunctions [28j. especially 
at small ri 2 ; (ii) by setting w A (r 12 ) = erf (Ar 12 )/r 12 (the 
"erf" adiabatic connection), the KS correlation energies 
from Eq. (|40|l have an error which is less than 4 mH for 
nuclear charges Z > 2; (iii) again with the "erf" adia- 
batic connection, we found that when the reference sys- 
tem corresponds to some A = [i between zero and A p h ys 
[as in Eq. (|21() ] our correlation energies are one order of 
magnitude better for ji > 1 jr s . 

The correlation potential of Eq. originally pro- 

posed by Overhauser [33| , has been also used to solve the 
effective Eqs. ll-iOl for the uniform electron gas (UEG), 
yielding to a very accurate de scrip tion of the short-range 
part of f{r 12) at all densities [23. A more sophisticated 
effective potential, based on a self-consistent Hartree ap- 
proximation, extended such accuracy to the long-range 
part of the UEG f{r 12) at metallic densities [24j. Other 
simple approximations for w^a (t*i 2 ) in the UEG have also 
been proposed and tested |25|. 

VII. CONCLUSIONS AND PERSPECTIVES 

We have presented the ideas concerning a theory 
based on the spherically and system-averaged pair den- 



sity f(r 12), and we have suggested to combine it with 
DFT to obtain system-adapted correlation energy func- 
tional. So far, the method has been tested for the He 
series [E| and for the uniform electron gas [23L I24I [25| , 
yielding promising results. 



In order to completely develop the approach presented 
here, many steps are still to be performed and many is- 
sues are to be addressed. Among them, the most rele- 
vant ones concern the construction of better approxima- 
tions for the effective electron-electron interaction that 
enters the formalism, and the implementation of a self- 
consistent scheme to combine the Kohn-Sham equations 
with the correlation energy functional arising from our 
approach. Last but not least, with the approximations 
tested so far our approach works very well for the short- 
range part of f(r 12), so that the combination with mul- 
tideterminantal DFT pH l34j (in which only the short- 
range correlations are treated within DFT) is also very 
promising and deserves further investigation. 
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